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We consider an electron constrained to move on a surface with revolution symmetry in 
the presence of a constant magnetic field B parallel to the surface axis. Depending on 
B and the surface geometry the transverse part of the spectrum typically exhibits many 
crossings which change to avoided crossings if a weak symmetry breaking interaction 
is introduced. We study the effect of such perturbations on the quantum propagation. 
This problem admits a natural reformulation to which tools from molecular dynamics 
can be applied. In turn, this leads to the study of a perturbation theory for the time 
"^h dependent Born-Oppenheimer approximation. 

a 



1 Introduction 

Recent advances in experimental physics have made it possible to produce two dimensional 
conducting surfaces of mesoscopic size. In such devices, the mean free path often exceeds 
the system size so the electron motion is ballistic and quantum coherence effects play a 
crucial role. This gives a motivation to strive for a complete understanding of the quantum 
mechanics of corresponding processes. In particular, conducting carbon "nanotubes" which 
are more or less uniform cylinders, belong to the family of surfaces that are nowadays ex- 
perimentally within reach ||. Since their discovery, lots of studies have been devoted to the 
elucidation of the spectral and transport properties of such devices, in a variety of situations 
and approximations - see, e.g. 0, ER, WFAS|| and references therein. Nanotubes of differ- 



ent types can be combined, and also coupled to other carbon structures such as fullerene 
molecules |[KBR|| , producing a variety of cylindrical surfaces . 



In this paper we study a quantum propagation in an "imperfect nanotube" subject to 
a constant magnetic field parallel to the tube axis within a simple model. Our model as- 
sumption is that a single electron is confined to a surface of revolution with slow variation 
of the radius along the revolution axis. Moreover, we assume that the rotational symmetry 
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is weakly violated, either by an impurity or by an external field. In other words, the used 
idealization amounts to neglecting the atomic structure of the tube as well as the interaction 
between the electrons, but taking into account the gross shape of the device. Our aim is 
to study the propagation of the electron along such an imperfect nanotube in the homoge- 
nous magnetic field by means of the time-dependent Schrodinger equation, starting with an 
initially localized wave packet, and paying a particular attention to the transitions between 
angular levels caused by the symmetry breaking perturbation. 

To understand the peculiarities of this quantum mechanical problem, it is useful to review 
briefly its classical counterpart; this is done in Section 2. The first question in the quantum 
case is the meaning of the fact that the electron is confined to a surface. The most natural 
approach, to our opinion, treats the surface as a limiting situation of a thin hard-wall layer. 
This idea goes back to [ |dC 1| , |dC2| , [To| and requires a renormalization in which the transverse 
contribution to the energy - blowing-up in the limit - is removed. One gets in this way an 
additional curvature-dependent term, in general attractive, to the potential. For the sake 
of completeness recall that there has been another recent work treating particle motion on 
revolution surfaces — see \ ~K\ , [Xr| , |1V1V| | and references therein. The last two papers aim at 
solvable models of compact surfaces (neglecting the curvature-dependent term), while 
treats the Schrodinger and wave equations on noncompact cylindrical surfaces without a 
magnetic field from the PDE point of view. 

Having thus found the Hamiltonian of our quantum system, we can analyze its spectral 
properties. When the rotational symmetry is preserved, we can perform (using a suitable 
gauge) the partial-wave decomposition. We can compute the angular part of the spectrum 
which depends on the actual cylinder radius varying along the tube axis. This brings to 
mind analogy with the molecular dynamics in which each angular state corresponds to an 
"electronic" level and the longitudinal coordinate measured at the axis corresponds to the 
one-dimensional "nuclear" configuration coordinate. Furthermore, when the rotational sym- 
metry is broken by a perturbation, the above analogy remains valid and we may invoke the 
time-dependent Born-Oppenheimer approximation to describe the propagation along cou- 
pled angular levels according to ||H1 , HJ1 , |HJ2| . Recall that the theory in molecular systems 
involves a small parameter which is given by the mass ratio between "electrons" and "nu- 
clei" . In our mesoscopic system, it is replaced by the parameter e defined as the inverse of 
the lenghtscale over which the variation of the radius of the nanotube takes place. Note, 



however, that we cannot directly apply the theory of [HI, HJ1, HJ2| in our perturbative 



context and a modification is needed as we shall explain below. 

The result of the analysis presented below gives a complete and rigorous description to 
the leading order of the wave function when the dynamics makes the electron go through a 
region where a perturbation couples angular levels. The basic picture is as follows. As long 
as the perturbed angular levels along the trajectory remain well isolated, the components of 
the wave function referring to the corresponding eigenstates are unchanged, to the leading 
order. When the unperturbed angular levels display a crossing or an avoided crossing, 
transitions between the perturbed angular eigenstates may become non-negligible as in the 



mentioned molecular analogy - see RHJ1|| . We are going to consider precisely the situation 
where the transition amplitudes are of order one, under perturbations of order -y/e. In such 
cases, an initial wave function having a nonzero component in a single angular eigenspace 
before the (avoided) crossing splits into the corresponding angular eigenstates according to 



2 



the Landau-Zener formula, again to the leading order. 

We have said already that despite being based on the paper [|HJ1|| , our description is 
not a direct application of the molecular time-dependent Born-Oppenheimer approximation. 
Indeed, in the rigorous derivations of this approximation, the "electronic" spectrum and the 
eigenstates, i.e. the spectrum at fixed coordinate along the rotation axis in our setting, are 
taken as given data, and the approximate solution to the molecular Schrodinger equation 
is constructed from this information - see [[Hl| , [HJ1] , |HJ2|1 . In our situation, by hypothesis, 
we only have access to that spectrum in a perturbative sense, and thus we need to develop 
a perturbative version of the time-dependent Born-Oppenheimer approximation that only 
requires knowledge of the leading terms of the Rayleigh-Schrodinger perturbation series. 
This is done in Section 5, where the main technical result of the paper is stated in Theorem 
|5T7| . We believe that it is of an independent interest. 

Precise statements of our results require a certain amount of notation and are therefore 
given below in Proposition |4.2j and Theorems . 



2 Classical Mechanics 

Let us start by describing the classical dynamics of the system. We consider a particle of 
mass m and charge e constrained to move on a smooth surface S with revolution symmetry 
around the axis OX in a homogeneous magnetic field B = Be x , B > 0, parallel to this axis. 

Using cylindrical coordinates, the surface is characterized by the smooth positive real 
valued function M 3 x h- » R(x) G such that 

x = x 

y = R(x)cos(6) } (2.1) 
z = R(x) sin(#) 

where (x, 6) E M. x S 1 . The squared length element on S is ds 2 = (1 + R' \x) 2 )dx 2 + R(x) 2 d0 2 , 
so the corresponding metric tensor g%j(x, 9) is given by 

= ( 1+ o W2 #(,))■ < 2 ' 2 > 

Using the circular gauge, we express the vector potential at the surface as 

A(r) = 1b A r = ( - sin(fl) | (2.3) 

This makes it possible to compute the Lagrangian function of the system 
L(r, r) = -mr 2 + erA(r) 

= \m [x 2 {\ + R\xf) + R 2 {x)6 2 ) + eBR ^ d . (2 . 4 ) 
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The system is integrable: we find that the momentum Pe = ^ and the kinetic energy T are 
two constants of motion 



_,,/ \ a eBR 2 (x) . r , 

p e = mR 2 (x)9 + ^-L , (2.5) 

T = lm(x 2 (l + R'(x) 2 )+R 2 (x)e 2 ) . (2.6) 



2 

Using ( |2.5| ) to express 9 as a function of x in Q2.6p , we deduce 



T = \mx 2 {l + R'{xf) + ^(jfi^-peeB + ^R 2 (x) 

= \mx 2 {l + R\xf) +V(R(x)). (2.7) 
The effective potential 3 R \— > V^(i?) G M + admits a unique minimum at i?o such that 

R« = \It^ and II (2 ' 8) 




Note that if pg = 0, the potential V(R) is harmonic on M.+ . From these considerations 
we deduce easily, in particular, that particle motions in the simplest case characterized by 
x(t) = correspond either to (x(t),9(t)) = (xq,9o) for any initial conditions (xq,9o), or 
to (x(t),9(t)) = (xq,9 — —t), where =: uj c is the cyclotronic frequency, for any initial 
conditions (xq, 9 ), or finally to (x(t), 9(t)) = (xo, 9o + ut), where ui is any constant, for initial 
conditions (xq, 9 ) such that R\xq) = 0. In case that R'(xo) ^ 0, the first two motions are 
stable, whereas in the last one the stability depends on the local properties of R around xq. 
In a similar way one can treat the general case with x(t) ^ 0. The motion is governed by 
the effective potential determined by the shape of S, and the potential minima correspond 
to the points where the angular motion has the cyclotronic frequency. 

Furthermore, notice that the addition of a supplementary exterior potential W, depending 
on x only, does not affect the functional dependence of po and its value remains independent 
of time. It is just the second constant of motion which is changed at that in the sense that 
the total energy E = T + W is now constant. 

Finally, let us also give the corresponding Hamiltonian function of the system for a future 
purpose. With p x = || we compute 

In the sequel we shall consider our charged particle to be an electron, e = — \e\ < 0, and use 
the rational units in which lei = m — 1 as well as h — c — 1. 



3 Quantum Mechanics 

Consider now the same system within quantum mechanics. For the purpose of this section, 
the function R : R — > M + defining the surface S is supposed to strictly positive and C 3 
smooth; later we shall impose stronger requirements. 
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The state Hilbert space of such a system is thus L 2 (S). To construct the Hamiltonian, 
however, it is not sufficient to replace the classical variables in ( |2.9|) by the corresponding 
canonical operators. The most natural way of quantization consists of taking a particle 
confined within a cylindrical layer built over S and squeezing its thickness to zero — cf. 
||dCl| , |dC2| , [Tc| |. One has to renormalize the energy in the limit, of course, subtracting the 



blowing-up part corresponding to the transverse motion. 

In the absence of the magnetic field, one arrives in this way to the Hamiltonian which 
equals to a sum of the respective Laplace-Beltrami operator (times 1/2 in our units) and the 
curvature-induced potential V(x) — — ~ (f?i(x) _1 — Q2(x)~ 1 ), where Qj(x), j = 1,2, are the 
principal curvature radii at the given point. The second part is of a purely quantum nature 
and has no classical counterpart. In the present case the locally elliptical intersection of S 
with the normal plane has the radius Qi(x) = R(x) , while for the intersection with the axial 
plane we find 

/ \ (i + R'jxyy/ 2 

= ^ — ; (3-1) 

the signs of Qi, q 2 coincide if both the osculation radii point the same side of the surface. 
Consequently, the curvature-induced potential equals 

V(x) = _ * ( t + ■ (3.2) 



8R{x) 2 V (1 + R'{x) 2 fl 2 



To express the kinetic (Laplace-Beltrami) part, — ~ |gr| 1 ^ 2 di\g\ 1 ^ 2 g^dj, we use ( j2.2|) and the 
corresponding contravariant tensor on S, 

(gV(x)) = Cd + y (3.3) 



R(x) 

The Hamiltonian in the presence of the magnetic field is then obtained by replacing the 
angular momentum operator pg = —idg by pg — A(x)R(x) where A(x) := Ag(r); it acts as 

rr 1 » R(x) a 1 ( , BR(x) 2 ^ 2 

H = = o x — = o x H ~, — — idg 

2R(x)^l + R'{x) 2 y/l + R'{x) 2 2R(x 



R{x)R"{x) 



8R{x) 2 V (1 + R'(x) 2 ) 3 / 2 
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1 + TT^T-bKTH (3-4) 



on an appropriate domain in L 2 (M x S 1 , R(x) \J 1 + R'(x) 2 dx d9). Due to the rotational sym- 
metry it has a simple partial-wave decomposition; its H m component is obtained replacing 
— idg by its eigenvalue m. In this way the spectral analysis of H is reduced to a family of 
one-dimensional Sturm-Liouville problems. Also the magnetic term has a natural meaning: 
we have 

, , , , , BR(x) 2 <&(x) , , , , ^ 

A(x)R(x) = — ±-L = ^ = cj>(x) , (3.5) 

where <fi is the magnetic flux value measured in the standard units (27r) _1 , or the number 
of flux quanta passing through the cross section of the cylinder. 
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It may be convenient to get rid of the weight factor replacing by an operator H on 
L 2 (IR) ®L 2 (5' 1 ). This is achieved by the unitary transformation ip i— > R 1 I 2 {1 + R' 2 ) 1 ^^. The 
only term in (|3.4j) which changes at that is the first one: by a straightforward computation 
we find 



H = - d x — . — XON d x + — xn H — 



-j +V 21 (x) + V 22 (x) (3.6) 



2(l + i?'(x) 2 ) a 2 J R(x) 2 

with 

^w=-^ri + „^» i <3.T) 



8R(x) 2 V (1 + R'{x) 2 f/ 2 
and 

rr / \ f ^ 2 _ 7 R' 2 R" 2 R" + R{R'R"' + R" 2 ) \ 

V22[X) ~ V 8R 2 (1 + R' 2 ) 8 (1 + R' 2 f + 4fl(l + R' 2 ) 2 ) [X) ' 1 ' 

Spectral properties of the Hamiltonian are influenced by the geometry of S. Suppose, 
e.g., that the latter has asymptotically constant radius, limui^co R{x) = Rq . In the absence 
of the magnetic field the problem is similar to that of a locally deformed Dirichlet strip 
BURS] , |KV| | (it is simpler, however, unless a mode-coupling perturbation is introduced). In 



the s-wave part the effective potential V 2 \ creates a potential well when S is locally squeezed 
and a barrier in case of a protrusion. For higher partial waves and non-zero magnetic field, 
of course, the effective potential consists of several competing contributions. 



4 Quantum Propagation 

Our main interest in this paper is not so much the spectrum of the Hamiltonian ( [3.4|) but 
rather the way in which an electron propagates over the surface of the cylinder. We will be 
particularly interested in the limiting situation when the radius modulation is gentle. This 
is conventionally described by means of the scaling transformation x i— > ex considering the 
asymptotic behaviour as e — » 0. This can be considered as a semiclassical limit since e — > 
means that the wave packet size becomes ultimately much smaller than the length scale of 
the radius variation. 

It is clear from the preceding section that the effective potential Vi = V21 + V22 is then 
dominated by the first term. Moreover, the operators ( |3.4j ) and ( |3.6| ) coincide in the leading 
term, which will be in the following the object of the investigation. We write its action as 

H(t) = 4 d 'TT7vM a - + ^ + uk) + (4 ' 1) 

on a suitable domain of L 2 (R) ® L 2 (S' 1 ) where R(x), V\{x) = R'(x) 2 are smooth on K and 
e) is smooth on 1 x [— e ,eo], for some eo > 0. Introducing an i?-dependent operator 
h(R) for R e R; by 

hiR) -M-' d ' +E r-)' < 4 - 2 ' 
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on a suitable domain of L 2 (S' 1 ), we can regard H(e) as an operator on L 2 (M, L 2 ^ 1 )) that 
we write as 

H(e) = -^ d * 1 + e l Vl (x) dx + ^ e) + H{R{X)) - (43) 
The spectral analysis of h(R) is straightforward and yields a family of simple eigenvalues, 

1 / BR 2 ^ 2 



nihill)) {A„(/?).,;eS}: { — [„ + — ) .neZ}. (J... I) 

with the corresponding eigenvectors 

<p n {6) = exp{in6)/V2^ , n G Z. (4.5) 

Note that the eigenvalues A n (i?) correspond to the classical effective potential V(i?) in fl2~~ 
with n G Z in place of pg. For n/mwe have 

A^)-A m (^) = ^(^ + £) (4.6) 



i? 2 



so that 



\ n (R) = \ m (R) <£> n + m<0 and R = R n , m =\l-±— i. (4.7) 



B 



Moreover, 



URn,m)= m ) 2 >0 . (4.8) 

Hence any pair of levels (A n (i?), X m (R)) with n + m < exhibits one and only one crossing 
as R varies, whereas other pairs never cross. The crossing points are well separated, 

{/?,,,„ : (n.m) 6 M". n m < 0} = { \j ^ ■ k G N* } . ( 1.9) 



with w -g = i? n _(fc +n ), n G N, and the values of the different pairs of levels crossing at 



for k fixed, are also well separated since 
wn 5(2n + fc) 

A„^ ni _ (fc+n) j - — . (4--WJ 

We note also that \ n (R) — A_ n (i2) = Sn. 

Thus, depending on our choice of function i?(x), the spectrum of h(R(x)) may display 
real or avoided crossings of an arbitrary width. Our aim is to adapt the techniques developed 
in | HJ1|| to describe the propagation of Gaussian wave packets (in the variable x) through 



these (avoided) crossings and, in particular, the splitting of the solution among the different 
angular levels X n (R(x)) involved. In particular, we want to allow an e dependent definition 
of the shape of our tube; it will then turn out that the natural scale for the phenomena we 
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want to describe is 8 = y/e. We henceforth adopt 8 as our small parameter and consider 
smooth functions R(x,S) defined on M x [— 8o,8q\. This means, in particular, that both the 
function V\ and the operator h will depend on both x and 8 in a smooth fashion. 

However, the above described model can exhibit no transitions because of the rotational 
invariance due to which passages between different levels A n are forbidden. To get a nontriv- 
ial result, we perturb therefore our system by introducing a real valued potential 8W(x, 9, 8), 
which is smooth on R x S 1 x [Sq, 8q] and violates the symmetry. For example, we can add 
a constant electric field in the direction d = sm(a)e z + cos(a)e x , where a (jL 1m. As a 
consequence, we lose integrability of the system on the classical level, whereas in the quan- 
tum setting transitions between the different perturbed eigenstates become possible. By 
assumption, when considered as a (bounded) operator on L 2 (S' 1 ) for (x,8) fixed, the oper- 
ator 8W(x, 9, 8) does not commute with h(R(x,8)), and therefore it perturbs the spectrum 
a(h(R(x, 8)). For the time being, let us keep the general form 8W(x, 9, 8) for the perturba- 
tion and describe the differences and similarities of the present case in comparison with the 
paper 



HJ1 



We introduce the operator g on (a suitable domain of) L 2 (R, L^S 1 )) by 

g(x, 8) = h(R(x, 8)) + V 2 {x, 8) + 8W(x, 9, 8) (4.11) 
so that the perturbed full Hamiltonian reads (with a slight abuse of notation) 

H M=-l S + 6>V 1 (x,S) a ' + 9(X ' S) - (412) 

Without loss of generality, we can assume that J sl W(x, 9, 8)d9 = by modifying V^x, 8) 
if necessary. We require the different potentials introduced so far to be smooth so that the 
following regularity hypothesis is fulfilled 

HO: The operator g is strongly C°° in (x,8) inlx [— 8o,8q\. 

We want to approximate the solutions to the Schrodinger equation in a suitable time scale, 

i8 2 ^ = h(8)4>, (4.13) 

for t in a finite time interval, as 8 — ► 0, for initial conditions of a "coherent state" type, 
which we shall describe in detail below. 

The first difference in comparison with |HJ1| comes from the fact that the kinetic term 



gives rise to a perturbed Laplacian 

8\ 1 _ 8\ 2 8\ VfaS) 

2° x l + 8*V 1 (x,8) x ~ 2° x+ 2° x l + 8±V 1 (x,8) x 



S -^d 2 x + R(x,d x ,8) (4.14) 



where 



R ^ = - S ° TTmks-y 2 " (-^ 1 + ^(^) 1 (-^*).(4.15) 
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We assume 



HI: 

sup V} k \x,6) <oo, fc = 0,l. (4.16) 

x6M,|5|<<5 

The factor S 8 in front of the operator R makes it possible to show that the influence of this 
term is negligible on the propagation of Gaussian states, so that the approximation given in 
[HJl l remains valid. This claim is the main result of this section and will be made precise 



in Proposition 4.2 below. 



The second difference in comparison with ||HJ1|] is that unless we have and explicitly 



solvable situation - and such are rare - we do not know in general the exact eigenvalues 
and eigenstates of the operator g(x,5). However, the approximation derived in |[HJ1|1 is 
constructed on the basis of this exact knowledge. A way out is to use an incomplete in- 
formation coming from the perturbation theory. Our second result, Theorem |5j], stated in 
section 5 says that it is enough to know the first few terms in the perturbation series in order 
to construct an approximation that describes the propagation, even in presence of avoided 
crossings, and that the result is as good as the one derived in [HJ1| . 



The rest of this section is organized as follows. We proceed with the description of the 
ingredients needed for our approximation, in analogy with |[HJ1|| , assuming that we know 



the exact diagonal form of g(x,5). Then prove that the modification R of the Laplacian 
does not affect the validity of this approximation. The next section will be devoted to the 
perturbative aspects mentioned above. 

We will denote by (j, n (x, 5) the eigenvalue of g(x, 5) such that fi n (x, 8)—X n (R(x, 5)) — > as 
<5 — > 0, for x such that R(x, S) far from R n ^ m . The corresponding eigenvector will be denoted 
by $ n (x, 5). If R(x, S) lies in a neighborhood of R ntm , we will denote by Ha(x, 5) > Hb{.%i 8) 
the almost degenerate perturbed eigenvalues with corresponding eigenvectors $a(x, 5) and 
5). The reason for such a convention is that the unperturbed eigenvalues X n (R(x, S) 
may or may not cross, are therefore the labeling of the /i's in terms of the indices n and m is 
not straightforward. Let Q n (x, 5) be the one-dimensional spectral projection of g(x, 5) corre- 
sponding to yU„(x, 5) in the first case and P(x, 6) be the two-dimensional spectral projection 
of g(x, 5) corresponding to Ha(x, 5) > ^b{ x i $) i n t ne second case. 

The situation we will study is that of avoided crossings of minimum width of order 8. 
Without loss of generality, we can assume the avoided crossing to occur around x = 0. More 
precisely we suppose that: 

H2: The eigenvalues /^(x, 5) and hq(x,5) are such that (/i^ — ^bY~^{0} = (0,0) in a 
neighborhood of (0,0) and mf xe i(^(x,8) — /ig(x,<5)) = c\5\ > for 5^0, where c is a 
constant and I is a small interval containing 0. 

We also set 

g l} (x,5) = g(x,5)P(x,5) (4.17) 
g ± (x,5) = g(x,5)(I-P(x,8)). (4.18) 



We know from [|H2| that locally around (0, 0) there exists an orthonormal basis, denoted 
as {ipx(x, S), ij)2{x, 5)}, of P(x, 8)L 2 (S 1 ), which is regular in (x,S) around (0,0). It is con- 
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structed in the standard Gram-Schmidt way: we choose an orthonormal basis {V'ljV^} of 
P(0,0)L 2 (S 1 ) and set 

" POT <« 9 > 

, r ^ (I-|0i(x^))(0 1 (x^)|)P(x,^ 2 

W ' } " ||(I-|0 1 (x,5))(0 1 (x,5)|)P(x,^ 2 ||- l4 - /Uj 
There exists a (x, 5) independent unitary transform U such that in the orthonormal basis 

if> j (x,S) = U<f> j (x,S), j = 1,2, (4.21) 
the matrix g\\ (x, 5) takes the form 



9\\(x,S) = gi{x,S) + V(x,S) 

P(x, 5) j(x, 5) + ia(x, 5) 

"f(x,S) — ia(x,S) —(3{x,5) 



+ V{x,5) (4.22) 



where V(x, 5) = trace(g(x, 5)P(x, 5)) is a regular function of (x, 5) around the origin and 

P(x,6) = hx + b 2 5 + 0(2), (4.23) 
7 (x,e) = c 2 5 + 0{2), 
a(x,5) = 0(2), 
V(x,8) = O(0), 

where b\ > 0, c 2 > 0, b 2 G K, and the following shorthand is used: 

0( m ) = O ((x 2 + 5 2 ) m/2 ) . (4.24) 

In order to get rid of the (^-dependence in the leading order of /3(x,5) in (|4.22|) , we 
introduce new variables, 

x' = b 1 x + b 2 8, 5' = c 2 5, t' = bl/c 2 2 t. (4.25) 
In terms of these variables, the Schrodinger equation ([4.13 ) for 

<j)(x', t') = ip(x(x', S'),t(t')) (4.26) 
becomes 

i5' 2 ^<P(x',t') = -^ r - ^— ^0(^,0 + C ^g{x{x',5'),5{5'mx',t') (4.27) 

in the limit 8' — *> 0, with 

V;(x',8')=V 1 (x(x',8'),8(6'))/4, (4.28) 

9ll (x(x',S%6(6')) = ( X ' 5 } f x ,\+0(2) + V{x(x',5'),5{5')), (4.29) 
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where V(x(x r , 5'), 5(5')) and V{(x',6') are regular in (x',5') around (0,0) and 0{2) refers to 
x' and 5'. We introduce the fixed parameter r = c\jb\ > and henceforth drop the primes 
on the new variables. We assume that gi(x, 5) has the form ( |4.22| ) with the following local 
behavior around x = and 5 = 0: 



P(x,6) 
7(x,<5) 
S(x,S) 
V(x,5) 



rx + 0(2) , 
r5 + 0(2), 
0(2), 
O(0) 



(4.30) 



with r > 0. 

Let us next describe the building blocks of our Born-Oppenheimer states. 

We begin with the definition of the semiclassical "nuclear" wave packets denoted as 
ifj(A, B, h, a, T], x). It comes from |[H3j| ; we have specify it for our setting where x G R. 
Since |[H3|| provides a detailed discussion of these wave packets, we refrain from proving all 
their properties here. 

We assume a G R, i] G R and h~ = 5 2 > 0. Let us stress that while the last symbol is 



useful when adapting the results of |[H3|| , it has nothing to do with the Planck's constant. 
We also assume that A and B are non-zero complex numbers that satisfy 



ReA B = 1 . 

This condition guarantees that Re BA~ X is positive, since ( Re BA 



(4.31) 



-i \- x 



\A\ 



Our definition of tfj(A, B, H, a, 77, x) is based on the following raising operator 



A(A,B,h,a,riY 



2h 



B(x 



A{-ih d --ri) 



(4.32) 



Definition: For the index j = 0, we define the normalized complex Gaussian wave packet 
(modulo the sign of the square root) by 



(p (A, B, h, a, rj, x) = 7r" 1/4 /T 1/4 A' 1/2 

x exp { - B A' 1 (x - a) 2 /(2h) + irj (x - a) /h) . 

Then for any positive integer j we define 

1 



(fj(A, B, h, a, 77, •) 



; A(A, B, h, a, n)*) J (p (A , B, h, a, 77, • ). 



(4.33) 



(4.34) 



Remarks: 1. For A = B = 1, H = 1, and a = 77 = 0, the <fij{A, B, h, a, 77, • ) are just the 
standard harmonic-oscillator eigenstates with energies j + 1/2. 

2. For each A, B, h, a, and 77, the set { ipj(A, B, H, a, 77, • ) } is an orthonormal basis for 
L 2 ' 



3. The position and momentum uncertainties of the <fj(A, B, h~, a, 77, ■ ) are *J~Q + l/2)h~ \A\ 
and a/ (j + l/2)h \B\, respectively. 



11 



4. When we solve approximately the Schrddinger equation, the choice of the sign of the 
square root in the definition of ifo(A, B, h, a, rj, • ) is determined by continuity in t after an 
arbitrary initial choice. 

5. Defining the scaled Fourier transform to be 

= (2irh)- 1/2 [ ^/(x)e-^ x/h dx, (4.35) 

J R 

then 

[Tn<fi(A,B,h,a,r,,-))(0 = (-i) 1 e~ iria/h Vi (B, A, h, r], -a, £). (4.36) 
We also define 



Vb ( x , 5) = V{x, 5) ± ^P 2 (x, S) + 7 2 (x, 5) + a 2 (x, 5) , (4.37) 
where x G IR, 5 > 0. Let a c (t) and r] c (t) be the solutions of the classical equations of motion 

j f a c (t) = V c (t), (4.38) 

±rf(t) = -VV c (a c (t),6), C = A, B , (4.39) 
at 

±S c (t) = V c (tr/2-V c (a c (t),5), 

with initial conditions 

a c (0) = 0, (4.40) 
r/ c (0) = V °(5), 

where 

rj°(5) = r] + O{5), rf > 0, (4.41) 
S c (0) = 0. 

The error term 0(6) depends here on whether C is A or B. In case of isolated eigenvalue /i n , 
V c = V n = /v 

We further introduce complex numbers which are defined by means of classical quantities. 
Let A c (t) and B c {t) be the solutions of the linear system 

jA c {t) = zB c (t) (4.42) 
±-B c {t) = i V c{2 \a c (t),6)A c (t) 



where a c (t) is the solution of ( |4.38|) and ( f4.40|) , with initial conditions 



Ay) = a , 
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It is easy to see that these quantities actually describe the linearized classical flow around 
the trajectory (a c (t) , rf (i) . The above convention regarding C applies if fi n is isolated in the 
spectrum. The asymptotics of these classical quantities for small t and S are described in 
detail in Section 2 of [ |HJ1| . 



The determination of the "electronic" part of the Born-Oppenheimer wave packet (B- 
O states, for short) is as follows. Although the "electronic" Hamiltonian is independent of 
time, it is convenient, since we deal with the time dependent Schrodinger equation, to choose 
specific time dependent "electronic" eigenvectors. Since they may become singular when the 
corresponding eigenvalues are degenerate, or almost degenerate, we shall define them for t in 
the outer regime, that is when a(t) is far enough from 0. This outer regime is characterized 
by times t such that (see ||HJ1|| ) 

S 1 ^ < \t\ < T, £ < 1/3. (4.44) 

We shall have two sets of eigenvectors, denoted by &q(x, t, 5), where the label ± refers to 
positive and negative times. Of course, this distinction is irrelevant if we consider an isolated 
eigenvalue /i. 

Let r] c (t) be the momentum solution of the classical equations of motion( 4.38| ) and ( 4.4Q ) . 
The normalized eigenvectors $^(x, t,<5) are the solutions of 

(§t(x,t,S)\ (d/dt + V e (t)d x )^(x,t,5)) = (4.45) 

for C = A, B and t^O. Since the eigenvalues 5) and /ig(x, 5) are non-degenerate for 

any time t small enough, such vectors exist, are unique up to an overall time independent 
phase factors, and are eigenvectors of gi(x,5) associated with Ec(x, 5) for any time. More 
precisely, we define the angles (p(x, S) and 9(x, S) by 

P(x,5) = W 2 {x, 5) + 7 2 (a;, 5) + a 2 {x, 5) cos{6{x, 5)) (4.46) 
7(z,<f) = y/P 2 (x, 6) + 7 2 (x, 5) + a 2 (x, 5) sin(0(x, 5)) cos(<p(x, 5)) (4.47) 
a(x,8) = ^(3 2 (x, 5) + 7 2 (x, 5) + a 2 (x, 5) sm(6(x, 5)) sm(<p(x, 5)). (4.48) 

and construct static eigenvectors. Let 

&a(x,8) = e i,pix ' 5) cos(6(x,5)/2)^ l (x,5) + sm(6(x,5)/2)^ 2 (x,5) (4.49) 
$b(x,5) = e- i ^ x ' s) cos(6(x,5)/2)ilj 2 (x,8) - sin(9(x,5)/2)ip 1 (x,5) (4.50) 

be the eigenvectors of gi(x, 5) associated with Hc{x, 5), C = A, B, for n/2 < 9(x, 5) < n, and 

$+(x,8) = cos(6(x,5)/2)ij 1 (x } 5)+e- i,fiix ' 6) sm(6(x,5)/2)ilj 2 (x,5) (4.51) 
$+(s,5) = cos(e(x,5)/2)ilj 2 (x,5)-e i ' pi:c ' s) sin(e(x,5)/2)ilj 1 (x,5) (4.52) 

be the eigenvectors of gi(x, 5) for < 9(x, 5) < ir/2. The solutions of ([4.45 ) are of the form 
*${x,t,S) = ^(x,S)e iX c(^), [l>l (4.53) 

where \f(x, t, 8) is a real valued function satisfying the equation 

z^A±(x, t, 5) + i V c (t)d x \±{x, t, 5) + ($ c (x, 5) \rf{t)d x $ e {x, 5)) = . (4.54) 
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We can get an expression for and its derivatives as follows. We fix values of the indices 
and drop them in the notation. We introduce the new variable 

uj = x - a(t) (4.55) 

and the notation 

X r (u,t,5) = X(uj + a(t),t,5) , (4.56) 

$ r (u,t,6)=$(u + a(t),6). (4.57) 



In terms of these new variables, equation Q4.54|) for X r reads 

d 



d 

i—X r (uj, t,8) = - ( <& r (a>, t, 5) 



m $ r (u,t,6)) (4.5* 



with 



^ r (oo, t, 5) = v(t)d x $(uj + a(t), 5). (4.59) 
By integration we get 

Xr(u,t,S) = J 7]{t')(<S>{u + a{t'),5)\d x <5>{uj + a{t'),5))dt' + X r0 {uj,5), (4.60) 

where we are free to set the integration constant X r0 (u, 5) = 0. 

The "nuclear" wave function is localized around the classical trajectory in the semiclas- 
sical regime. In view of the genericity condition rf > 0, in the outer temporal region major 
part of the "nuclear" wave function will be supported away from the neighborhood where the 
levels almost cross. Hence we can introduce a cutoff function which does not significantly al- 
ter the solution and forces the support of the wave function to be away of this neighborhood. 
Let F be a C°° cutoff function, 

F:R+^M, (4.61) 

such that 

F(r) = l ... 0<r<l 

F(r) = ... r > 2 { > 

The wave functions we construct below in the outer regime will be multiplied by the regu- 
larizing factor 

FiWx-a'm/S 1 - 6 '), (4.63) 

where < 5' < f , for C = A, B. 

Remark: On the support of F the relation 

x = V °(5)t + 0(5 1 - 6 ' + t 2 ) (4.64) 
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holds true, and since r]°(5) = rf + 0(5), where rj > 0, we find that 
\x\ > c\t\, 



uniformly in 5. 



A Born-Oppenheimer state ift^fe, 5, t) is defined by 



(4.65) 



(4.66) 



F{\\x-a c {t)\\/5 l - 5 ') Vj {A c {t), B c (t), 5\ a c (t), ?? c (t), x )e isC ^ s2 ^(x, t, 5). 

It is a good approximation to the solution of the Schrodinger equation (|4.13|) as 5 — > 0, when 
R is absent and far enough of the crossing region, i.e. in the outer time regime (|4.44 ), as 
shown in ||HJ1|| . Proposition fO] below shows this is still true when R is present. 
In the inner time regime, characterized by the inequality (see [|HJ1|| ) 

- 6 1 -* < t < S 1 ^ , ^ < 1/3 , (4.67) 

we look for an approximation constructed by means of the classical quantities associated 
with the potential V(x, S), the average of fi (x, S) and (i B (x,5). Let a(t) and S(t) be the 
corresponding classical quantities satisfying the initial conditions 

a(0) = 0, 

r)(0) = rf , (4.68) 
5(0) = 0. 

It is suitable to use the rescaled variables 

V = (*-<W S (4.69) 

It is shown in ||HJi| that a good approximation ipi of solutions to ( [4.13|) in that regime, when 
R is absent, is given by 



My,s,S) = F(\\y\\6 5 ')exp [i 



S(5s) .rj{5s)y 



+ i 



S 



with 



X(y, s, 5) = f (y, s)^i(a(<fa) + 6y, 6) + g (y, s)tp 2 (a(Ss) 
where fo,go are complex-valued functions solutions to 



Sy,S) 



9 ( fa(y,- 



ds \ga(y,s) 



rfs + y 1 \ / f {y,s) 

1 -{rfs + y) ) \g (y,s) 



(4.70) 
(4.71) 
(4.72) 



The general solution to this equation is 

/ (l-D / 



fa(y,s) 

9o(y,s) 



CM 



(i-O 

2 



D ir 



27 ? U 



{-l + i)^{rfs + y) 



D 



I i)J^(rfs-y) 



(4.73) 



+ C 2 (y) 



D 



2n u 



-(l + i)/^( v °s + y) 



(l-M) 



+\D_ 



!2>j u 



"(I 



^(v°s + y) 
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The coefficients C\{y) and Ci(y) have to be determined by matching with the incoming 
solutions of B-0 type at the border of the inner and outer time regimes. 

In particular, assuming for defmiteness that the incoming B-0 state ipoi is associated 
with the index / for the "nuclear" component and the B level with the initial momentum 
rf{Q) = rj°, i.e. that 

tPoiix, t, 6) = (x, t,6), -T<t< -8 1 -* , (4.74) 
we have 

Ci(y)=0 (4.75) 

and 

C 2 (y)= - S-^wiAo, B -irA /r)°, l,0,0,y)e"^ exp (jL(y>-2y) 

x exp (i S ° ^ ~^ + -^(1 + 31n(2?7 ) + mr - A\nS)] , (4.76) 



5 2 Ar]° 

where Sq(S, — ) is real and can be computed explicitly - see |[H Jl|] . 



The analysis of ||HJ 1|| shows that in this situation, we get an outgoing solution given 



by a linear combination of B-0 states, with explicit coefficients, associated with the same 
index / for the "nuclear" component but with both levels. The initial momentum is chosen 
as 77 (0) = i]° — 2rS/rj° for the A level and the outgoing solution ipoo is of the form 

r(i + ^) 

provided <t <T, where 

X(5) = Tr/4 + S£(-,5)/5 2 + ^--r (l + 31n(277°) + lnr-41n<fi . (4.78) 

2t] u 



Here again, Sq-(5, — ) is real and can be computed explicitly from |[HJ1 



Moreover, the function obtained by pasting the approximations constructed in the outer 
and inner temporal regions is an approximate solution to the Schrodinger equation when the 
perturbation R of the Laplacian is absent (see [4.14] ). Similar explicit formulae are valid if 



the ingoing state is associated with the A level. Hence, the propagation through avoided 
crossings can be iterated. 

We are going to show that the perturbation of the Laplacian in ( (4.14j ) does not affect the 
propagation of B-0 states. The general strategy is simple: we write 

H(5) = -jA x + g(x,8)+R(x,d x ,5) 

= H (5) + R(x,d x ,5) (4.79) 

and denote by ^i(x, t, 8) the approximation given by ipoi^i^oo i n their respective time 
domains constructed in ||HJ1|| : 

ipoi(x,t,S) ... 

y l (x,t,8) = { tpi(x,t,5) ... -($!-€< t < (J 1 -^ (4.8O) 
ipoo(x,t,5) ... 5^<t<T 
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We define (i by 

= iS 2 dt%(x,t,S)-H(S)%(x,t,S) 

= iS 2 d t %(x, t, 5) - H Q (6)%(x, t, 5) - R(x, d x , 5)%(x, t, 5) 

= $(x,t,5) + $(x,t,5), (4.81) 

where is the error term controlled in ||HJ1|| by means of the following abstract lemma. 



Lemma 4.1 Suppose H(h) is a family of self-adjoint operators labelled by h > 0. Suppose 

that ijj(t, h) belongs to the domain of H(h), is continuously differentiable in t, and solves 
approximately the Schrddinger equation in the sense that 

ih^{t,h) = H(h)i>{t,h) + ((t,h), (4.82) 

where £(t, K), satisfies 

|| c(t, K) || < fi(t, h) . (4.83) 

Then 



e 



itHW/n^fy _ ^h)\\ < fir 1 / fx(s,h)ds (4.84) 



o 



holds true for t > and the analogous statement is valid for t < 0. 
Using the same lemma to estimate the norm of (I, we get 



Proposition 4.2 Under the hypotheses (H0)-(H2), the function ^fi(x,t,8) defined by ( \j-8Q) 
is for any T > an approximation to the solution if)(x,t,8) of the Schrddinger equation 
fcTTSj ) such that 

ij(x,t } 5) = %(x,t,S) + 0(5 P ) (4.85) 
holds in the L 2 (IR) sense for some p > and all t G [-T, T\. 
The proof of this technical proposition is given in the appendix. 

5 Propagation of Perturbed B-O States 

Let us now turn to the second indicated step and replace the above B-0 approximation 
by a construction making use of a perturbative knowledge of the exact eigenvectors and 
eigenvalues of the operator g\\ defined by fl4.22|) . In particular, this needs to be done for the 
quantities appearing in Q4.38|) , (|4.42|) determined by means of a classical potential given by 



an approximation of the spectrum of g\\ . We will show that it is enough to know the second 
order perturbation expansion in order to construct an approximation of the solution that is 
a perturbed version of our initial B-0 states and still describes accurately the transitions 
between the "electronic" levels. 
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In order to make some explicit formulae simpler and to stress the effect of the perturba- 
tion, we will assume in this section that both the operators h and W are 5 independent, i.e., 
we shall consider 

g(x, 6) = h(R(x)) + V 2 {x) + SW(x, 9) , (5.1) 

where V2 commutes with h whereas W doesn't. This means that h(x) is assumed to have 
a degeneracy at x = in the considered part of its spectrum that is lifted by W to the 
leading order in 5. This is the generic situation we set out to investigate when the avoided 
crossing results from a weak symmetry breaking violation of a true eigenvalue crossing. Note, 
however, that we are able to accommodate the general situation considered so far, at the 
cost of more complicated perturbation formulae. 

Let us state a simple lemma which is at the basis of our constructions and which says 
that an approximation of an approximate solution is an approximate solution. 

Lemma 5.1 Let H(5) be for all 5 G (0, 80) a self adjoint operator densely defined in a 
Hilbert space H, and let ip a {t,5) G TL, <p a (t,8) G H be time dependent vectors with the 
following property: there exist c,pi,p 2 > such that the relations 

|| e -iW^ oM - VaM)ll < c5 Pl (5.2) 

and 

\\<Pa{t,S)-Mt,8)\\ < cSP2 ( 5 - 3 ) 

hold for all t from an interval Ici and < 5 < 5q. Then 

e- iHms2 cp a {0, 5) - <p a (t, 5) < 3c5 min(pi ' P2) , 

e -^(«5)*/^ a ( 0) S) - <p a {t, 5) < 3c5 min{pi ' P2) . (5.4) 
Proof uses just unitarity of the evolution group and the Cauchy-Schwarz inequality. □ 

Our approximate B-0 states will require classical quantities defined by means of an 
approximation V c of the potential V c used in ( [L38Q , ( 4.42 ). We have to estimate the error 
induced by this approximation. In order to do that, we make use of Gronwall's lemma (see 
e.g. 0) that we recall below. 

Lemma 5.2 Let E be a Banach space, U C E be open, I be an interval of IR and f G 
C 1 (J x U;E) be such that there exists K > O with sytQ^^ixu HAa/^j ^)IU(b) < K. Let 
g : / x U — > E be continuous and such that there exists G > with 

sup \\g(t,x)\\ <G. (5.5) 

(t,x)eixu 

If a and f3 be C 1 maps from J — > U ( where J C I) satisfying for t G J 

a\t) = f(t,a(t)), (5.6) 

(3'(t) = f(t,(l(t)) + eg(t,P(t)), (5.7) 

then 

\\a{t) - /3{t)\\ < ||a(t ) -/9(to)||e x| *-* 01 +eG(e^-*°l - l)/K. (5.8) 
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For convenience let us recall here our definition ( |4.66| ) of a Born-Oppenheimer state 
ijjj(x,5,t) in the exterior regime: 

tf(x,8,t) = (5.9) 
F(\\x-a c (t)\\/5 1 - s ')^(A C (t), B c (t), 5\ a c (t), V c (t), x <S>±( Xl t, 5) . 

We want to compare ip^(x, 5, t) with an altered but similar definition based on approximate 
quantum and classical quantities for ipj(x, 5, t): 

^ C (x,S,t)= (5.10) 
F{\\x-7f{t)\\/5 l - 5 ')^{A c {t), B c (t), 6\a c (t), f(t), x )e rsC ^ s2 ^(x, t, 6) . 

All "tilded" classical quantities are generated by equations ( |4.38|, |4.42| ) with an approximate 
potential V c (x, S) in place of V c (x, 5). The vector 

$±( X) t, 5) = $±(x, 5)e a c(^) (5.11) 

depends on the approximate classical quantities through the phase and on an approximate 
normalized quantum eigenstate (x, 5). Note that we keep the same Gaussian function cpj 
to construct the "nuclear" wave packet. 

Our next goal is to apply Lemma |5.1| to estimate the errors in terms of the difference 
between V and V c . 

Lemma 5.3 The following inequality holds in the outer time regime for the L 2 (M) norm: 

\\^ ± (x,S,t)-^ C± (x,5,t)\\ < (5.12) 
c (\A(t) - A(t)\ + \B(t) - B(t)\ + \a{t) - a(t)\/5 2 + \rj(t) - V (t)\/5 2 

+^ sup {\rj(s)- V (s)\ + \V(a(s))-V(a(s))\ + sup \d x V (x)\\a(s) - a{s)\) 
° se[o,t] xe[5(«),o(«)] 

+ sup \F(\\x - H(t)||/5 1 -' 5 ')$±(a;, t, 5) - F(\\x - a{t)\\/ d 1 ' 5 ')^, t,6)\. 

(x,t,8) 

with some constant c. 

Proof: The index C being fixed in this context, it will now be omitted. Other irrelevant 
parameters will also be dropped in the arguments. Note that since the function F is smooth, 
we can write 

F{\\x - a^W/S 1 - 6 ') = F(\\x - a(f)||/<f 1 -*') + 0{{a{t) - a(t))/5 l - s ') (5.13) 
and that the L 2 (S' 1 )-norm of the vectors $^(x, t,<5) equals one. Since 

V(a) - V{a) = V(a) - V(a) + V(a) - V(a) (5.14) 
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and 77 and rj are uniformly bounded, we infer 

S(t) = J (rf(s)/2 - V(a(s))ds^ = 3{t)+ (5.15) 

0\\t\sup l\rf(s)-rj(s)\ + \V(a(s))-V(a(s))\+ sup \d x V (x)\\a(s) - a(s)\ 

V se[0,t] y ase[5(s),a(s) 

Then we compute 

(pi(A,B,d~ 2 ,a,rj,x) - ipi(A,B,5 2 ,a,r],x) (5.16) 
= e^"^/* 3 5, 5 2 , a, 0, re) - B, 5 2 , a, 0, x) 

+(pt(A, B, 5 2 , a, 0, x) L^ir^/P) - e i{v-v)(x-°<)/s 2 \ jvi 



From Lemma 3.1 in ||HJ1|| we learn that as A — > A and 5 — > B 

(pi(A,B,5 2 ,a,0,x) = <fi(A,B,S 2 ,a,0,x) 

+ 0(\A-A\ + \B-B\ + \a-a\/5) 



holds in the L 2 (M) sense, which takes care of the first term. Then we note that the L 2 norm 
of the remaining term is equal to 

||^(A B, S 2 , a, 0, x )e l ^ a -~ a ^ - cp^A, B, S 2 , a, V -?], x) || (5.17) 
= \\ Vl {B, A, 5 2 , 0, -a, ^e^-^/- 52 - e^^^^B, A, 5 2 , rj-% -a, a;) || 
= 0{(a-a)/5 2 + (v-v)/S 2 ) 

by using Plancherel formula, the properties of the (pj under Fourier transform, \\<Pj\\ = 1 and 
the above lemma again. Then, gathering these estimates and using the facts that A(t) and 
B(t) are uniformly bounded, we get the result. □ 

In order to use the just proved lemma, we see that it is necessary to approximate V c to 
an error of order o(5 2 ) and to show that this induces errors of the same order in the classical 
trajectory (a c , fj ) and errors of order o(l) in the linearized classical flow (A c , B c ). Moreover, 
the corresponding eigenstates should be at most at a distance o(l) from $ c . 

When we consider times away of the matching regime, i.e. r < \t\ < T, where r is 
independent of 5, it is easy to show the following result, just by using Gronwall's lemma and 
regular perturbation theory. We thus omit the proof. 



Lemma 5.4 Let the time interval (r, T) be such that the solutions to (U-38j ), M.4-0Q satisfy 
the condition 

i {a c (t) \r<t<T, < S < 5 } = P , (5.18) 
where the corresponding potential 

V c {x,S) =/i C (x,5), xeP, (5.19) 
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is the nondegenerate eigenvalue of g(x,5) corresponding to ^c{x,S). Let 

V c {x, S) = //g(x) + 6f4(x) + 6 2 fj£(x), xe P, (5.20) 

be the second-order perturbation expansion for /x c (x, S) . We define a c , fj c , A c , B c , S c as above 
with the conditions 

a c {r) = a c (r) + o{5 2 ) , f(r) = rf{r) + o{5 2 ) , (5.21) 

A c (r) = A c (t) + o(l) , B c (t) = B c (r) + o(l) , (5.22) 

5 c (t) = S c (t)+o(5 2 ), (5.23) 

and 

$ c (x,t,5) = ^ c (x,0)e fXcix ' t ' s) , (5.24) 

where Xc(x,t,8) is given by ( U-54) with ^c(x,0) in place of^c(x,5) an d 

\e{x,T,8) = X c {x,T,5) + o{l). (5.25) 



Then there exists a solution ifj(x,t,5) to the equation ( \4-l$ ) such that 

ip{x,t,5) (5.26) 
= F{\\x-lf{t)\\/^ 5 ')^M C {t), B c (t), 5 2 , a c (t), rf(t), x )e i ~ sC ^ 5 ^c(x, t, 5) + o(l) 

holds true in the L 2 -sense and for all t < t < T . 

Remark: We have the familiar explicit formulae 

f4(x) = fi c (x,0), (5.27) 
l4(x) = ($e(x,0)\W(x)<S> c (x,0)) , (5.28) 
t4(x) = -($c(x,0)|^(x)(/i(x)- /U c (x,0)) r - 1 ^(x)$c(x,0)), (5.29) 

where the reduced resolvent is given by 

(M,)-/(,o,)^E ( ^'M . ,5,0) 

The above result has to be modified for times close to the matching regime, since in 
that case degenerate perturbation theory is required to define the potential. Indeed, the 
approximate potential chosen in the lemma diverg — > 0, so that Gronwall's lemma 
cannot be used as it stands. Let us find the modified potential from the perturbation theory. 

The two eigenvalues of g(x, 6) which are of interest to us, /^(x, 8) and /iz?(x, 5), are given 
by the spectrum of P(x, S)(h(x) + 5W(x) + V^x)) This operator is represented in the smooth 
orthonormal eigenbasis ( |4.21| ) by the matrix ( |4.22| ) , which we can expand to second order 



in 5 for any i in a neighborhood of the origin, since the projection P(x, 5) entering the 
definition of the basis ( f4.21| ) is regular. Hence we can write 

( x . ( P(x,$) j(x,8) +ia(x,6)\ ~ (K , 

9 ^={y(x,6)-i*(x,6) -P(x,S) ) +V ^ > ( 5 - 31 ) 
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where 



P(x, 5) = /3q(x) + 5(3 1 (x) + 5 2 (3 2 (x) + 0(5 3 ) = B 3 (x, 5) + 0(5 3 ) 

j(x,S) = j (x) + S ll (x) + S 2 l2 (x) + O(S 3 ) = G 3 (x,S) + O(S 3 ) 

<j(x,5) = a (x) + 5a 1 (x) + 5 2 a 2 (x) + O(5 3 ) = S 3 (x,5) + 0(5 3 ) 

V{x, 5) = V {x) + 5V 1 (x) + 5 2 V 2 (x) + 0{5 3 ) = V 3 {x, S) + 0(S 3 ) 

with the error 0(5 3 ) being C°° in x, and (see ( ]4.30| )) 

Po{x) = rx + 0{x 2 ) , Pi(x) = 0(x), 

l0 (x) = O(x 2 ), 7l (x) = r + 0(x) , 

a {x) = O(x 2 ), a l {x) = 0{x). 

Let us set 



(5.32) 
(5.33) 
(5.34) 
(5.35) 

(5.36) 
(5.37) 
(5.38) 



s(x, 5) = y/(B 3 (x, 5)) 2 + (G 3 {x, S)) 2 + (S 3 (x, 5))* (5.39) 
and define our (explicit) modified potential by 

V c (x, 5) = ±s(x, 5) + V 3 (x, 5) (5.40) 
where the sign is chosen according to the value of C. It is easy to check that by construction, 

V c (x, 5) - V c {x, 5) = 0{5 3 ) (5.41) 

as x — ► 0. As above, we employ tilde to mark the values generated by the modified potential. 
We only consider the dynamics for positive times, the other case being similar. 

To define the perturbed classical trajectory, we will start integrating Newton's equations 
from a positive t (5) = 5 K , for some 2/3 < k < 1, using as initial condition the explicit 
asymptotic expansion given in Corollary 2.1 of | |HJ1| |: 

Corollary 5.5 In the outer regime 5 — > 0, t — > 0, |t|/d> — ► oo and t 3 /5 2 — > ; we have 



t 1 



a»(t)= - d x V 3 (0,5)- + V °(S)t±- 



+ 



2 

5 2 \n\t\ 

2(#)f 



5 2 



St 



S 2 In S 

4(^(5)) 2Vi ' y " !>> 2(r)°(S)) 2 



;i + 21n(2ry°( ( 5))) 



+ 0(t 3 ) + 0(5 4 /t 2 ) 

The asymptotics for rf{t) in the same regime is obtained by termwise differentiation of the 



above formulae up to errors 0(t 2 ) 

The choice of t (5) ensures that 

a(t ) = a(t ) + o(5 2 ) , 
Wo) = v(to)+o(S). 



0(5 A /t 3 ). 



(5.42) 
(5.43) 



Whereas the error is small enough for the position, it is not the case for the momentum. 
Hence we resort to energy conservation in order to determine the momentum with a sufficient 
accuracy. 
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Let us first note that due to the uniform boundedness of the force induced by the poten- 
tials V c and V c , there exist a r small but independent of S and constants < C\ < C 2 < oo, 
such that as long as t G [— r, r], 

Ci<7j c (t)<C 2 , (5.44) 

and similarly for rj . 

The unperturbed energy is given by 

E C (5) = (f(t)) 2 /2 + V c (a c (t), 5) = ( V C (5)) 2 /2 + V c (0, 5) (5.45) 
= E c (5) + V c (0,5)-V c (0,5) 

where the perturbed energy E C (S) is explicit. Hence 



r] c (t) = ^2(E c (S)-V c (a c (t),S)) > (5.46) 
holds for t G [— r, r], and we can define rj c (t) by energy conservation so that 



rf(t) := p{EC{5)-VC{7iC{t),5)) (5.47) 
= V c (t) + O (V c (a c (t), 5) - V c (a c (t), 5)) + O (v c (0, 5) - V c (0, 6) 

= r] c (t) + (sup \V c (x, 5) -V c (x,5)\ +a c (t)-a c (t) 



Thus using formula ( |5.47[ ), we make an error in r] c of the same order as the error we make 
in a and V c . 

Next we turn to the approximations A c (t) and B c (t). They are defined as solutions to 
the system ( |4.42j ) with V c in place of V c and initial conditions at t — ±to, given by 



AB(t)\_f A(0) 



B b (t) J yB^Tsign^rAiO)/^)) 



(5.4J 



It remains finally to consider the perturbed eigenvectors $c{%,t,5) (where we dropped the 
index referring to the sign of t). The restrictions to the support of F mentioned in lemma 
and the estimate (|5.44j) imply that if we impose the condition 



1 -<*'-«> , (5.49) 
we can write 

x = a(t) + 0(<5 1 - <5 ') > ct{l + OiS 1 ' 5 '^) > cS K (5.50) 

for some positive constant c, and the same estimate is true with a replaced by a. 

Hence in the considered regime the eigenvalues fi c (x, 0) of P(x, 0)(h(x) + ^(x)) display 
a gap that is at least of order x = 0(5 K ) - see the behaviour ( f4.30| ) - and we call the 
corresponding eigenvectors Xc( x )- We define our perturbed static eigenvectors by 

®c(x,6) = X c(x) (5.51) 
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and similarly, the phase corresponding to time dependent perturbed eigenvectors ^c(x,t,8) 
- in view (|5.11|) - by 



(5.52) 



X Cr (uj,t,5) = i / dsrj(s)(<&c(u + a.(s),5)\d x ® c (u + a(s),6)) 

'MS) 



where we used the new variables flOp and ( fL6Cp . 

The next lemma tells us that our definitions of (a c (t),rj c (t)), (A (t), B c (t)) and 3>c{x, S, t) 
are accurate enough for our purpose. The proof can be found in appendix. 

Lemma 5.6 With the definitions above, there exists a positive r such that for allt e [to(S), r] 
we have 



a c {t) = a c {t) + o{5 2 ), 

f{t) = V c (t) + o(5 2 ), 

A c (t) = A c (t) + o(l), 

_ B c {t) = B c (t) + o(l), 

$(x,t,5) = $(x,t,6) + o(l). 



(5.53) 
(5.54) 
(5.55) 
(5.56) 
(5.57) 



Hence, with the definitions made above, we have a perturbed B-0 state given by ( |5.10| ) 
that is explicitly expressed by means of perturbation theory in 5 (modulo finding the solution 
of the classical equations of motion, of course) and which yields an approximation of the 
solution to the Schrodinger equation Q4.13 ) for finite time intervals as 5 — > 0. In particular, 
putting together our results, we get the following statement. 



Theorem 5.7 Adopt the hypotheses (HI) and (H2) and assume the behaviors (\j.3(\) . Sup- 



-c± 



pose that 2/3 < k < 1 and r is as in the above lemma. Let ipj (x,8,t) with \t\ > 5 K be a 



perturbed B-0 states according to (\5.1Q) constructed by means of the approximate quantities 
considered in lemma \5.4\ if \t\ > r and in lemma \5.(\ if 5 K < \t\ < t, subject to the condi- 
tion that all classical quantities agree at the instants t = ±r. Let ip(x,5,t) be a solution to 



equation fy.lSj ) with ip(x, 5, —T) 



ip(x,5,t) =ij)f (x,6,t)+o(l), 
holds as 5 — > for all —T < t < —5 K , while 

ip(x,5,t) = -e-* r ^°4>f + (x,5,t) + e-^ /4 "° 



x,5, — T). Then 



{5.5i 



iir 




p i\{5) 



ipj (x,6,t) + o(l) (5.59) 



holds for all 5 K < t < T , with X(5) given by ( |^. 78 ). 



Remark: 

It is possible also to give an explicit approximation of the wave function in the inner time 
regime, —5 K < t < 5 K , in terms of quantities coming from perturbation theory. However, 
this temporal region being so short, it is not crucial for most applications to have a detailed 
approximation there. 
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6 Appendix 



Proof of Proposition |4.2| : It is enough to show that the norm of in (|4.81|) is small and 
apply to Lemma P . The expression ( |4.15| ) together with HI show that we only need to 
control the effect of p = —i5 2 d x and p 2 = (—i5 2 d x ) 2 on ^i(x, t, 5), since for any tp we have 



\\R{x,d x ,5)ij;\\<C{5*\\ P 2 i)\\+6< 



(6.1) 



First consider the outer temporal region and the form ( 4.66 ). We know from the computa- 
tions in [ |H3|| that 



-T))(fi{A, B, 5 2 , a, r/, 
- v) 2( Pi(A B, 5 2 , a, T], 
Moreover, we estimate 



\B\5y/l + 1/2. 



\B\ 2 5 2 v / (6l 2 + 6l + 3)/4. 



pF (\\x-a c (t)\\/5 1 - s '^ 
p 2 F(\\x-a c {t)\\/5'- 5 ') 



< c 2 6 2 ^ 



(6.2) 
(6.3) 

(6.4) 
(6.5) 



where the constants Ci, c 2 depend on F only. Away from the crossing region, the "electronic" 
eigenvectors are well defined and smooth in (x,5). Hence we only need to consider what is 
going on in the neighborhood of x = to get an upper bound on the effect of p and p 2 on 
the eigenvectors Q^(x,t,5) given by (|4.53| ). We drop the indices and consider 

${x,t,5) = e iA(x '*' 5) $(x,<5) (6.6) 

where $(x, S) denote some static eigenvectors and X(x, 8, t) the corresponding real valued 
function defined by fl4.54|) . We compute 



(6.7) 
(6.8) 



d x $(x,t,6) ■ 

d 2 x $(x, t, 5) = e iX( - x ^ [d 2 x <$>(x, 5) + 2(id x X(x, t, S))d x Q(x, 5) 
{id 2 x \(x, t, 5) - (d x X(x, t, 5)) 2 }$(x, 5)] . 

As rj, ipj, d x ipj, dfyj, d x ipj are all 0(0) as (x, t) — > (0, 0) in the support of F, we have 

d x $(x,6) = 0(d x 9(x) + d x <p(x)) + O(0) , (6.9) 
%$(x,8) = O ((d x 6(x)) 2 + (d xV (x)) 2 + d 2 J(x) + d 2 xV (x)) +O(0) (6.10) 

in the norm of the "electronic" Hilbert space. In the expression ( f4.60p for A, we first check 
by inspection that in all cases 



($\d x $) = O(d x ip) + O(0) 



(6.11) 



(see, e.g., (3.50) in | [HJ 1| | ) since all functions of 6 and are uniformly bounded and, moreover, 
the factor of d x ip is a function of 9 only. Hence, by a further differentiation we get 



d x (<s>\d x ®) = o (dl<p + d x <p + d x <pdj + d x e) + o(o) , 
&Md x $) = o(dlv + d 2 x vd x e + d x vd 2 x e + (d xV ) 2 + d 2 x v 
d x( pd x e + d 2 x e + (d x ef + d x e + d x y) + 0(0) 



(6.12) 
(6.13) 
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It remains to estimate d x ip and d x 9. We have 

<p(x,5) = arcta,n(cr(x,5)/j(x,5)) (6.14) 
provided j(x, 5) is different from zero. Hence using ( |4.30| ) we get 

^ 5) = >y*(x,6) + v*(x,6) ' (6 - 15) 

so that with the help of estimates of the type \^y\/ a/7 2 + cr 2 < 1 we arrive at 

(616) 

By similar operations we eventually obtain 



V a/7 2 + ^ 2 / V 7 2 + o- 2 / 



and 



9 3 = p [ + ) ( dlldxQ + gg}%7 + <9 2 7<9x7 + d 2 x ad x g 
" Va/t^T^/ V 7 2 + ^ 2 

c ( {d xl ?d x a + {d x afd xl + {d xl f + (d x af \ 

Assuming further that 

\\x\\=0(5 K ), £<2/3</c<l-f<l, (6.19) 
we get from the behaviour (|4.30| ) in this region 

= O (^) , (6.20) 

1 

6 

Then we consider 

6(x, 5) = arccos I ^ 6 2 = ) . (6.23) 

W0 2 (x,5)+i 2 (x,5) + a*(x,5)J 



%P = O ( \ 1 , (6.21) 
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By computing derivatives and estimating as above, we easily get 



dj 

die 



o 



o 



d x (3 



(3 2 + 7 2 + a 2 
d 2 J 



+ o 



+ o 



%l + d x a 



+ O 
+ O 



/? 2 + 7 2 + <x 2 , 

d x p(d x r/ + d x a) 
{(3 2 + 1 2 + a 2 )^{ 1 2 + a 2 ) 



v /(/5 2 + 7 2 + a 2 )( 7 2 + ( j 2 )^ 

v /(/3 2 + 7 2 + a 2 )( 7 2 + ( 7 2 ) i 

{d x l) 2 + (d x a) 2 + d x ~jd x a 
V / (/9 2 + 7 2 + ^ 2 )(7 2 + ^ 2 ) 



(6.24) 



+ o 



V/5 2 + 7 2 + a 2 ^ 
Using dQ5p , $ K C > |*| > 5^, and (|6~T9| ), we thus find 



d x 6 



die 



o 
o 



$2(1-0 

1 



(6.25) 

(6.26) 
(6.27) 



$3(1-5) 

Gathering the different pieces, we obtain for the derivatives of A in the regime just described 

X(x,t,S) = 0{t/5 1 - K ) = 0{l/5 1 - 2R ) , (6.28) 
d x X(x,t,S) = O (l/5 3 - 2 ^- 2K ) , (6.29) 
d 2 x X(x,t,6) = (l/5 A -^' K ) . (6.30) 

so that we obtain the following estimates for the derivatives of the vector $(x, t, S) 

d x ${x,t,6) = 0(l/6 2 - 2 ^) , (6.31) 
%$(x,t,6) = (l/5 4 -^) . (6.32) 

We are now in a position to estimate the effect of p and p 2 on the B-0 states in the outer 
time regime: 

= \\{pF)^ e + F(p Vl )& + F Vl (p<b c )\\ 
< c(5 1+s ' + \\(p- V c ) Vl \\ + \ V c \+5 2 \d x <5> c 



< c(l)(6 1+s ' +5B e + \ V e \ + 5 2 ^ . 



(6.33) 



We have already used above the fact that \r] c \ is uniformly bounded as 5 and t go to zero, 
and the same is true for B c - see Lemma 2.1 and Proposition 2.2 in ||HJ1|| . Finally we get 
in the outer temporal regime 



||Wfll<c(0 



(6.34) 



as 5 — > 0, where c(l) is some constant independent of time. By similar manipulations we 
also get in the same regime 

lb 2 Vf || < c(/) . (6.35) 
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1 



Note that the non-vanishing term comes only from the action of p on the Gaussian state 
(pi, which yields essentially rj c as expected, whereas the contribution from the derivatives of 
the "electronic" eigenvectors and cutoff function vanish. From the definition of R we get a 
supplementary 5 4 which more than compensates for the denominator 5 2 appearing in ( 4.84 ) 

/ \\R{x,d x ,S)7pf{x,t,6)\\dt<c{l)6 A . (6.36) 

We now need to perform the same type of analysis on the approximate wavefunction 
ipi(y, s, 5) given by ( (4.70[ ) adopted in the inner temporal region. There we use the variables 
( |4.69| ) so that the relations 

d y = 5d x and p = —i5d y (6.37) 

have to be employed to compute the derivatives of the different pieces in the definition of 
ipi(y, s, 5). In this case we need to show that 

\\R(x,d x ,5)My(^t),s(t),8)\\dt 



1 



1 

6 



\R(6y + a(s5),6dy,5)if} I (y, s,6)\ 2 5dy \ ds 



5-5 



1 /2 

< t^7 sup ( [\R(6y + a(s5),6d y ,5)My,s,5)\ 2 5dy\ (6.38) 

is of order S p for some positive p as 5 — > 0. As above, we denoted at that the norm in the 
"electronic" Hilbert space by a modulus. The estimates ( |6.4| ), fl6.5| ) remain valid and we 
have 

\ p 2 e iv(ss)y/s\ = \rf( s 5)\ < C, (6.40) 

since rj(t) is uniformly bounded in the inner temporal regime. Noting that x = a(5s) + 5y = 
0(5 l ~^), we also get from the regularity of the orthonormal basis {ipi(x, 5),i/j2(x, 5)} around 
(0, 0) that 

\p^(a(5s) + Sy,S)\ = 0(5 2 ) , (6.41) 
b^H^ + ^l = 0(5 4 ) (6.42) 

for j = 1,2. Finally, the functions fo(y,s) and go(y,s) determined in (|4.73|) to (|4.76| ) and 
their derivatives can be estimated using the following remark. Up to phases, these functions 
are given as products of a Gaussian, a polynomial in y, a parabolic cylinder function times, 
and a factor 1/S 1 ^ 2 coming from the normalization of the function (pi. Asymptotically, these 
parabolic cylinder functions, their first and second derivatives are of order 0((s + \\y\\)°), 
0((s + \\y\\)) and 0((s + \\y\\) 2 ), respectively, where s = 0{5~^). Hence we can write 

|/o(lM)| < Pi(y)e-y 2 ^ 2 6^ 2 , (6.43) 
\pfo(y,s)\ < P 2 (2/)e^ 2 /2|Ao| 2 r i/2+i-^ (644) 

\p 2 fo(y,s)\ < P 3 (y)e-y 2 ^ 2 5-^ 2 - 2 ^ (6.45) 



28 



where Aq is the initial condition ( |4.43| ) and Pj, j = 1, 2, 3, are polynomials in y the coefficients 
of which are independent of 5. They depend on I, the index of the chosen B-0 state. Similar 
estimates are valid for go in place of fo- Having 

p ^ = e iS(Ss)/S^ V (sS)y/S [(pF + i^)( /o ^ + g ^ 2 ) + F((p/ )Vl + (p<7o)V>2 

+/o(W>i)+ 00(^2))] (6-46) 

and the above estimates we can write 

< P 4 (l/)e- s,2/2|A,|2 r 1/2 (5 1+5 ' + l + 5 2 + 5 1 -^ (6.47) 

with another polynomial P4. Hence the right hand side of ( |6.38| ) can be further estimated 
to give 



1/2 

\ P My,s,6)\ 2 5dy\ <c(/)/5 1+ «. (6.48) 



-<5-e<s<<5-« 
By similar manipulations we also get 

1 

-<5-e<s<<5-« 



1/2 

f . 1+s sup { I \p 2 My,s,S)\ 2 6dy} <c{l)/5 1 ^. (6.49) 



We note that here the leading order contribution comes from the action of p on the phase 
e ir iv/ 5 which gives i]. The supplementary factor 5 A in (|6.1| ) yields the final estimate 



1 



\\R(x, d x , 6)Mvfa t), s(t), 5) \\dt < c(l)5 3 ^. (6.50) 



Hence the proposition holds with p = 3 — £. □ 

Proof of Lemma |5.6| : As noted above, we cannot directly use Gronwall's lemma as stated 
in the text. Hence we need to prove that the two evolutions stay close enough to each other 
between times to(S) and r, where r will be small but independent of 5 by a more refined 
analysis. We consider the index A and drop it in the notation. 

First, it is easy to check the following asymptotic properties as (x, 5) — > (0, 0): 

s(x, 5) - ^(3 2 {x, 5) + 7 2 (x, 5) + a 2 {x, 5) = 0{5 3 ) , (6.51) 

d x s(x, S) - ^v^M + fM+^M = 0( J— ) , (6.52) 

yx z + o z 

s(x, 5) = rVx 2 + 5 2 {1 + 0(x + 5)) , (6.53) 

d x s(x,S) = 0(l), (6.54) 

We collect some preliminary observations on the solution a(t) to the equation 

a(t)=-d x V(a(t),8) (6.56) 
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for t G [to, t] with initial condition satisfying ( |5.42| ). We can choose r > independent of 5, 
such that 

a(t) - a(t (6)) > c (t-t (5)), (6.57) 



for some c > and all t G [to(5), r\. This implies easily by means of (|5.42j ) that 

a(t)-a(t (5))> Cl (t-t (5)) (6.58) 

for all t G [to(5), t] with a uniform constant. Hence we can write 

x 2 + 5\ > c 3 (6 K + (t- t (5))) 2 , W t G [a(t), a(t)} . (6.59) 

Consider now the identities (dropping the 5 dependence in the arguments) 
5(t) - a(t) = d x V(a(t)) - d x V(a(t)) = 

d x y//3 2 + 7 2 + a 2 (a(t)) + d x V(a(t)) - d x s(a(t)) - d x V 3 (a(t)) = 

d x (^/W+l 2 ~T^(a(t)) - s(a(t))) + d x (V(a(t)) - V 3 (a(t))) ~ 

dls(9 t )(a(t) - a(*)) + d?V 3 (9 t )(a(t) - a(t)) , (6.60) 

where 6 t G (a(t),a(t)). Now the first order derivatives are of order <5 3 /(<5 K + (t — t (5)), 
whereas the second order ones are of order 5 3 / (5 K + (t — t )) 3 ~ see ( |6.55| ) and ( |6.59| ). 
Hence introducing d(t) = a(t) — a(t) we get an ODE of the form 

d'(t) = f(d(t),t)d(t)+g(d(t),t), (6.61) 
where we have the a priori bounds 

\f(d(s), s)\ds < c5 2 j 1/{5 K + (s- t )) 3 ds < c5 2(1 ~ K) , (6.62) 

to J tg 

and since we can assume without loss that 5 k + (t — t (5)) < 1 , 

/ W(s),s)\ < / — — — ^^<c5 3 (|ln(^)| + |ln(^ + (t-t )|) 

= 0{5 3 \n{5)). (6.63) 



Equation Q6.61|) is equivalent to 

d(t) = d{t ) + (t - t )d(t ) + [ ds [ du{f{d{u),u)d(u) + g(d(u),u)) . (6.64) 

J to J to 

Let us denote 

D(t) = sup \d(s)\ . (6.65) 
se[t ,t] 

We deduce from the above bounds 

\d(t)\ < \d(t )\ +(t- t )\d(t )\ +c [ ds D(s)6 2( ~ 1 -^ + cd 3 \ \n(8)\ < 

Jto 

c(\d(t )\ + \d(t )\ + 5 3 \\n(5)\+ [ ds D(s)5 2( - 1 -^) (6.66) 

Jto 

30 



and, as D is not decreasing, 

D(t) < c {\d{t )\ + \d{t )\+5 3 \ln{5)\ + f dsD(s)5 2{l - K) ) 

J to 

< c(\d(t )\ + \d(t )\ + 6 3 \ In(5)| + D(t)5 2 ^) . 
Since S 2 ^^ — > 0, we find that 

£>(*) <c(|d(t )| + |ci(t )|+5 3 |ln(5)|). 
Plugging this into ( |6.64| ) finally yields 

d(t) = d(t ) + (t- t Q )d(t ) + O(S 2 ^- K \\d{t )\ + \d(t )\) + 5 3 \ bi(S)\). 

As an immediate consequence of this result and (|5.47| ) we have for any t G [to{8), t] with our 
choice of to(^) and initial conditions ( |5.42| ) 

a{t) - a(t) = o(5 2 ) , (6.70) 
rj(t) - V (t) = o(5 2 ) . (6.71) 



(6.67) 
(6.68) 
(6.69) 



Turning to (A(t), B(t)) and their approximations, we first note that by [|HJ1| , p. 102] we 
have with our choice of t (S) 



(A(t ) 
\B(t ) 



A(t ) 
B(t ) 



oil). 



Then we consider the equation (equivalent to (|4.42|) and (f4.43|) ) 



(A(t)\ = (A(t ) 



+ 



to 



O i 
idlV{a{t)) 



Ms) 

Bis) 



(6.72) 



(6.73) 



\B(t)J \B(t Q) 

and a similar one for the approximations with the tilded symbols everywhere. Introducing 



A(t) 



/ A(t)\ _ ( A(t)\ 
\B{t)) \B{t)) 



we compute 



A(t) = A(t ) + 

' / 

id 2 x V{a{s)) - dlV{a{s)) 



(6.74) 



to 



A{s) 
B(s) 



ds + 



to 



i 
id 2 x V(a(s)) 



A(s) ds 



But 



\B(t)J 



0(1) by jHJll , /* d 2 V(a(s))ds = £>(5 2(1 ~ k) ) and similarly for the untilded 



quantities. Hence using the same type of manipulations as above, we deduce 

\\A(t)\\<c(5 2 ^ + \\A(t )\\). 
It follows that 



(6.75) 



(6.76) 
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holds for any t G [to(5), r). 

In order to deal with the "electronic" eigenvectors we consider the perturbation series for 
the resolvent (h(x) + V^x) + 5W(x) — z)~ l when the argument z runs through the circle of 
radius |// (x, 0) — /j, b (x, 0)|/4 centered at any of the eigenvalue /i C (x, 0). Integration on this 
circle yields the eigenprojector Qj(x, 5), j = A, B, and the estimates 



Qj(x, 5) = Qj(x, 0) + 0(SW(x)/\fj, A (x, 0) - /i B (x, 0)|) = Qj(x, 0) + 0(5/ x) , (6.77) 
d x Q 3 (x,5) = d x Q 3 (x,0) + O(5/\fi A (x,0) - fi B (x,0)\ 2 ) = 0) + 0(<5/x 2 ) .(6.78) 



This, in turn, yields the following estimates on the eigenvectors 3>j(x, 5) of the perturbed 
operator h(x) + 5W(x): 



Now we consider one eigenvector Xji x ) an d drop the index j. We note here that eq. (3.58) 



so that X r (u, t (5), 5) = 0(5 2k x ) — > with S. On the other hand, using the fact that x( x ) 
is smooth and that rj(t) is uniformly bounded on [io(5),T] we find 



^(x,5) = Xj(x) + 0(6/x), 
d s $j(x,6) = d xXj (x)+0(S/x 2 ). 



(6.79) 
(6.80) 



in |[HJ1|| shows that 



X r (co,t,S) = Oit/S 1 -") 



(6.81) 




Jt (8) 

x<($(w, s, 5) + 0(6 /(a(t) + u)))\ (d x $(u + a(s), 5) + 0(6/(a(s) + ^) 2 )))> 




o(l) + 0(51n((t + u)/(t + cu))) + O(5(l/(t (5) + u) - l/(t + u))) . 



(6.82) 



Having u = O^ 1 - 5 ') and (|5^9| ), the error terms above can be estimated by 



o(l) + C?(51n(5) + 5 1 - K ) 



(6.83) 



which goes to zero as 5 —> 0. It follows then that 



A(x,t,5) - \{x,t,6) = o(l) 



(6.84) 



and in turn we get 



(6.85) 



which concludes the proof. 



□ 
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